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Abstract 

We employ a functional approach to investigate the confinement problem in 
quenched Landau gauge QCD. We demonstrate analytically that a linear rising po- 
tential between massive quarks is generated by infrared singularities in the dressed 
quark-gluon vertex. The selfconsistent mechanism that generates these singular- 
ities is driven by the scalar Dirac amplitudes of the full vertex and the quark 
propagator. These can only be present when chiral symmetry is broken. We have 
thus uncovered a novel mechanism that directly links chiral symmetry breaking 
with confinement. 
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More than thirty years after the identification of non-Abelian gauge the- 
ory as the appropriate framework to describe the strong interactions we still 
lack a satisfactory understanding of the confinement phenomenon. Isolated 
particles with nonvanishing colour charges have not been observed in nature. 
This fact is supposed to be encoded in the infrared structure of QCD. In the 
quenched theory a gauge invariant signature of confinement is the area law 
of the Wilson loop at large distances. This behaviour corresponds to a linear 
rising potential between static colour charges in the fundamental represen- 
tation of the gauge group. This signature has been unambiguously verified 
in lattice gauge QCD, see e.g. ref. [TJ. However, the dynamical mechanism 
that generates the Wilson potential is still elusive. Since the underlying long 
range interaction is provided by gauge dependent objects, this mechanism 
may have a different appearance in different gauges. In this letter we present 
such a mechanism for covariant Landau gauge QCD. 

The idea of 'infrared slavery', i.e. the notion that infrared singularities 
generate confinement, dates back to the early seventies E]. Speculations 
about the infrared behaviour of the running coupling at that time distin- 
guish two cases: (A) the coupling develops an infrared fixed point and (B) 
the coupling diverges at small or zero momenta. The last possibility has been 
considered as the driving mechanism for the generation of infrared singular- 
ities and quark confinement. Today there is ample evidence from functional 
approaches that the first possibility is realised in the Yang-Mills sector of 
Landau gauge QCD: the running coupling freezes out at small momenta 
[U El [7] . However, this does not imply that the Green's functions of Yang- 
Mills theory are finite in the infrared: Instead, depending on the number of 
ghost and gluon legs, many of the one-particle irreducible Green's functions 
are indeed singular in the limit that all external momenta vanish [7] . In this 
letter we demonstrate that these singularities also induce a corresponding 
singularity in the dressed quark-gluon vertex. As a consequence the running 
coupling defined from this vertex diverges at vanishing renormalization scale, 
and a linear potential between massive quarks is generated thereby realizing 
infrared slavery via an infrared singular quark-gluon vertex. 

Besides confinement the other fundamental property of infrared QCD is 
dynamical chiral symmetry breaking, i. e. the nonperturbative generation of 
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quark masses from dynamical interaction with gluons. Lattice QCD tells us 
that the chiral phase transition and the deconfmement transition take place 
at a similar temperature [8]. This suggests a relation between confinement 
and dynamical chiral symmetry breaking, which, however, has not been clar- 
ified in detail yet. In this letter we present, based mostly on an analytical 
calculation, an explicit mechanism that links these two phenomena. 

The functional approach we employ for our investigation is given by the 
tower of Dyson-Schwinger equations for the one-particle irreducible (1PI) 
Green's functions of QCD [9j [TOj [TT] . A great advantage of this continuum 
based formulation as compared to lattice QCD is the analytical access to the 
infrared behaviour of these functions without finite- volume effects necessarily 
present in lattice calculations. As has been demonstrated for the propagators 
of three-dimensional Yang-Mills theory extremely large lattices are manda- 
tory to access the small momentum region where the infrared asymptotic 
scaling can be identified [12]. The computational costs for a corresponding 
investigation in four dimensions exceeds current possibilities. This problem is 
even more severe for higher n-point functions. This underlines the necessity 
for a continuum based approach to the Green's functions of QCD. 

To begin with we shortly summarise previous results for the infrared 
behaviour of the Green's function in Landau gauge Yang-Mills theory. Gauge 
fixing in this framework is performed by the standard Faddeev- Popov method 
supplemented by auxiliary conditions such that the generating functional 
consists of an integral over gauge field configurations that are contained in 
the first Gribov region. (The feasibility of this method has been justified in 
a framework employing stochastic quantisation [13] ■) The resulting Dyson- 
Schwinger equations for lPI-Green's functions have been solved analytically 
in the infrared to all orders in a skeleton expansion (i.e. a loop expansion 
using full propagators and vertices) [TJ. Here we are interested in the case 
where all external momenta go to zero. Choosing all momenta proportional 
to each other and requiring for the largest one p 2 <C Aq CD a self-consistent 
solution of the whole (untruncated) tower of DSEs is given by [7J 



where T n,m stands for the infrared leading dressing function of a Green's 




(1) 
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functions with 2n external ghost legs and m external gluon legs. The expo- 
nent k is known to be positive [2]. An important property of the infrared 
solution (TjQ) is the fact that it is generated by exactly those parts of the DSEs 
that involve ghost loops [7] . In other words: the Faddeev-Popov determinant 
dominates the infrared behaviour of non-Abelian Yang-Mills theories. Thus 
an infrared asymptotic theory can be obtained by 'quenching' the Yang-Mills 
action, i.e. setting exp(— Sym) — > 1 in the generating functional [13]. The 
solution of this asymptotic theory is given by the power laws ([1]). Interest- 
ingly, this limit is a continuum analogue of the strong coupling limit of lattice 
gauge theory. 

Examples of the general solution ([T]) are given by the inverse ghost and 
gluon dressing functions, r 1,0 (p 2 ) = G^ 1 (p 2 ) and T 0,2 (p 2 ) = Z _1 (p 2 ), respec- 
tively. They are related to the ghost and gluon propagators via 

The corresponding power laws in the infrared are 

G(p 2 ) ~ bT K , Z{p 2 ) ~ (p 2 ) 2 * . (3) 

Since k is positive one obtains an infrared diverging ghost propagator, a 
behaviour which is necessary to ensure a well defined, i. e. unbroken, global 
colour charge [15]. In Landau gauge an explicit value for k can be derived 
from the observation that the dressed ghost-gluon vertex becomes (almost) 
bare in the infrared [TBI E3 HBJ- One then obtains k = (93 - VT201)/98 « 
0.595 [5], El [19], which implies that the gluon propagator vanishes in the 
infrared. A direct consequence of this behaviour are positivity violations 
in the gluon propagator and therefore the confinement of transverse gluons 

[HED]. 

A further important consequence of the solution (pQ) is the qualitative 
universality of the running coupling in the infrared. Renormalisation group 
invariant couplings can be defined from either of the primitively divergent 
vertices of Yang-Mills-theory, i.e. from the ghost-gluon vertex (gh — gl), the 
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three-gluon vertex (3g) or the four-gluon vertex (4g) via 



* 9h - 9 \p 2 ) = ^G\v 2 )Z{ V 2 ) P ~° Cl /N c , 

« 3 V) = f^V)] 2 Z\p 2 ) p2 ~° c 2 /N c , 

« 4 V) = 9 ^\ V 2 )fZ\v 2 ) P -° c 3 /N c . (4) 

Employing the DSE-solution ([1]) it is easy to see that all three couplings 
approach a fixed point in the infrared. However, the explicit value of the 
fixed point, ci )2j 3 /N c , may be different for each coupling. For a bare ghost- 
gluon vertex one obtains a gh ~ gl (0) ~ 8.92/ N c [5J; the other couplings have 
not been determined yet. Equations P| underline that indeed the fixed 
point scenario (A) is realised in the Yang-Mills sector of QCD. In turn this 
explains the existence of the power law solutions (CD): the theory becomes 
approximately conformal in the far infrared. 

We now present a significant extension of this analysis to the quark sector 
of quenched QCD. To this end we first choose the masses of the valence 
quarks to be large, i.e. m> Aq CD . The remaining scales below A qCD are those 
of the external momenta of the Green's functions. Here we are primarily 
interested in the case where all external momenta go to zero. Then, without 
loss of generality, one can choose all external momenta to be proportional to 
one scale p 2 , which is small compared to all other scales in the theory, i.e. 
p 2 <C Aq CD . One can then employ Dyson-Schwinger equations to search for 
selfconsistent solutions in terms of powers of p 2 . 



integrals with an infrared power-law enhancement are themselves dominated by in- 
frared integration momenta. We illustrate this with the following simple example, where 
a > 1 is a real number. Consider 



Then, asymptotically for small p, 

1 1 



Hp) 



a - lp"" 1 

and it is the lower integration limit that contributes. Then we just need to subdivide the 
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Figure 1: Dyson-Schwinger equation for the quark-gluon vertex. In the first 
line we show the full equation [21]. The second line shows the lowest or- 
der of a skeleton expansion, that already reveals the correct power counting 
in the infrared. All internal lines in the diagrams represent fully dressed 
propagators. 

We apply this method to the Dyson-Schwinger equation for the full quark- 
gluon vertex, given diagrammatically in the first line of Fig. [U In the second 
line the higher n-point functions of this equation have been expanded to 
lowest order in a skeleton expansion in terms of full propagators and vertices. 
The dressed quark-gluon vertex T M can be decomposed in a basis of twelve 
tensor structures, which are given explicitly e.g. in [22]. Due to momentum 
conservation the vertex depends only on the two external quark momenta 
(p^, g M ) or three Lorentz invariants {p 2 ,q 2 ,p ■ q). To analyse the infrared 
limit of this vertex in the presence of only one external scale p 2 <C Aq CD , we 
can set g M = 2p^ without loss of generality. This leaves us with only four 
possible tensor structures which can be denoted by 

4 

r,(p)=^5>0° 2 ) G ; ( 5 ) 

i=l 

integration interval 

A Ai A 2 A„ 

/ = /+/+-/ 

Ai A„_i 

and it is always the lower (0) limit of the first subinterval that generates the power-law in 
p. Therefore, for low-p, the q integration is infrared dominated. 

Another way of putting it is that the integral is divergent with p = 0, and keeping finite p 
regulates it. Then the integral is dominated by momentum scales of order the regulator. 
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with £j being Lorentz and Dirac scalar functions and 

Gl = l», Gl = p,, Gj = j%, Gl=h», (6) 

where we have normalised the momentum, p^ = p^j ^fp 2 , to ease power 
counting. Note that G 1 ^ and G 3 ^ have an odd number of 7-matrices, whereas 
G 2 ^ and G^ have an even number. Therefore £2,4 7^ only if chiral symmetry 
is broken. 

The internal loops of the skeleton expansion in Fig. [1] are dominated by 
loop momenta similar to the external scale p 2 due to the denominators of 
the loop propagators. Thus for p 2 <C Aq CD the dressing functions of internal 
propagators and vertices can be approximated by the power laws given in 
fll]). We therefore use ([3]) for the scaling of the ghost and gluon propagators. 
The ghost-gluon vertex scales as a constant in the infrared [THl dZJ [TS] . The 
internal quark propagator lines can be written as 

_ qj + M(p 2 ) ij>Z f Z f 

S{P) ~ p 2 + M 2 (p 2 ) Zf{P ~M ' (?) 

for momenta p 2 < Aq CD , M = M(p 2 ->• 0) ~ A QCD and Z f = Z f (p 2 -> 
0). Hereby we assume that neither the quark mass function M(p 2 ) nor the 
wave function renormalisation Zf(p 2 ) is singular in the infrared [23]. This 
assumption will be justified below by selfconsistency arguments. From the 
two tensor structures in (J7J) the scalar piece turns out to be the leading one 
in the infrared. 

For the internal quark-gluon vertices we can employ the expression (131) 
with any internal momentum as argument: it contains all possible types 
of Dirac structures (vector, scalar, and tensor) and any more complicated 
dependence on external and internal momenta will generate the same powers 
of external momenta in ^..4 after integration (for dimensional reasons). To 
determine the infrared exponents of the quark-gluon vertex we employ the 
ansatz 

£i~(pY\ e 2 ~(pY 2 , £3~(pY 3 , u-( P 2 f 4 (8) 

for the scaling of the dressing functions of the quark-gluon vertex with mo- 
mentum. We then substitute this ansatz into the skeleton expansion of the 
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vertex-DSE and determine the exponents /3i..4 selfconsistently by matching 
the left and right hand sides. The resulting infrared solution is given by 



fo = -l/2-K, 



fafoe [(-1/2- k), (-k)] 
p A e[(-i/2-K),(i/2- K )]. 



(9) 



We have checked this analytic solution also in numerical calculations and 
found that the case 



is realised. The details of this calculation will be presented elsewhere [21]. 

There are two important remarks here: Firstly, similar to the DSEs in the 
Yang-Mills sector it is the diagram containing the ghost loop that dominates 
the right hand side of the equation in Fig. [TJ Thus also the infrared behaviour 
of the quark sector is dominated by effects generated from the Faddeev-Popov 
determinant. Secondly, the driving tensor structures of the solution (jHJ), (fTUj) 
are the scalar Dirac amplitude G 2 ^ of the quark-gluon vertex and the scalar 
(Zf/M) part of the quark propagator. Both structures are only present when 
chiral symmetry is broken. If they are absent, as in the case of restored chiral 
symmetry, one obtains a different solution. We analyse this case below. 

Before we discuss the solution (jUJ) further it is important to check that it 
persists to all orders in the skeleton expansion. Higher order terms in this 
expansion can be generated by inserting diagrammatical pieces like 



into lower order diagrams (all propagators are fully dressed). In ref. [7] we 
showed that insertions involving pure Yang-Mills propagators and vertices 
do not change the overall infrared divergence of a given diagram. This is 
also true for the three pieces involving quark lines: The first one introduces 
two quark-gluon vertices, two quark propagators and one gluon propagator 
plus a new integration into a given diagram. This amounts to an additional 
power of 



= -1/2 - K 



(10) 




(...) 



Z 2 (p2)-l-2«+2«-l+2 Z 2 



(11) 
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in the external momentum (p 2 ). A similar result is obtained for the other 
pieces. Thus diagrammatic insertions of quark lines do not change the 
infrared behaviour of a given diagram but introduce additional powers of 
Zf/M. This implies that the solution ([§]) is valid to all orders in the skeleton 
expansion. It is therefore also an infrared solution of the full vertex-DSE. 

We wish to add that the infrared divergence presented here has been 
found under the pretext that all external momenta of the vertex go to zero. 
However, a more detailed analysis of the vertex-DSE also shows, that the 
power law must be more general, (pl)~ 1 ' 2 ~ K in terms of the gluon momentum 
P3. The uniform divergence presented here is therefore a particular case since, 
once all momenta vanish, also the gluon momentum vanishes. The details of 
the more general analysis will be presented elsewhere [23] . 

An important application of our infrared power counting scheme concerns 
the four-quark 1PI Green's function H(p), which is given in Fig. [2] together 
with its skeleton expansion. From the terms of the skeleton expansion one 
obtains H(p) ~ {p 2 )~ 2 in the infrared. The well known relation 



between the static four-quark function H(p° = 0, p) and the quark potential 
V(r) therefore gives a linear rising potential by naive dimensional arguments. 
A more refined treatment, as described in [25J, leads to the same result. Note 
that already the first term in the skeleton expansion, i.e. nonperturbative 
one-gluon exchange displayed in Fig. [2J, generates this result. Since the fol- 
lowing terms in the expansion are equally enhanced in the infrared, the string 
tension will be built up by summing over an infinite number of diagrams. The 
coefficients of these diagrams will be investigated elsewhere [24J . 



Figure 2: The four-quark 1PI Green's function and the first terms of its 
skeleton expansion 




r 



(12) 
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A further interesting quantity is the running coupling a qg from the quark- 
gluon vertex. A nonperturbative and renormalisation group invariant defini- 
tion of this coupling is given by 

"«(P 2 ) = ^(i(p 2 )Zj(p 2 )Z(p 2 )~~, (13) 

where £1 dresses the 7 M -part of the vertex, see Eq. (jSJ). With the exponents 
P2 — Pi — —1/2 — k, as obtained in our numerical solutions, we find by 
power counting that the coupling is proportional to l/p 2 . Thus, contrary to 
the couplings (j3J) from the Yang-Mills vertices, we find this coupling to be 
singular in the infrared, i.e. 'infrared slavery' is realised. 

So far we found quark confinement in the chirally broken phase of quenched 
QCD. Our power counting scheme worked in this case, because only one ex- 
ternal momentum and no other small scales were present. This is also true 
in the chirally symmetric phase of (quenched) QCD with massless valence 
quarks. The solution (JTJ) for the pure Yang-Mills sector is still valid in this 
case. In the quark sector all tensor structures in the quark propagator and 
the quark-gluon vertex that violate chiral symmetry have to disappear, in 
particular £ 2 and £4 have to vanish identically. We are thus left with 

S(p) = ^Z f (p 2 ), T,(p)=zgY,UP 2 K (14) 

^ t=l,3 

for the quark propagator and for the quark-gluon vertex. On general grounds 
Zf(p 2 ) cannot be infrared singular and can be treated as a constant. We then 
obtain 

/3 1 = p 3 = -K (15) 

for the exponents in ([8]). This solution no longer leads to a confining poten- 
tial. One obtains H(p) ~ (p 2 ) -1 for the infrared behaviour of the four-quark 
function. This leads to a potential of Coulomb type 

V(r) = -^—r [ H(p° = 0, p)e ipr d 3 p ~ ^ . (16) 



(2^ 

Also the resulting running coupling from the quark-gluon vertex is no longer 
diverging but goes to a fixed point in the infrared 

a q9 (p 2 ) = ^e 1 (P 2 )zj(p 2 )z( P 2 ) ~ ^, (i7) 
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similar to the couplings from the Yang-Mills vertices. The restoration of 
chiral symmetry is therefore directly linked with the disappearance of infrared 
slavery and confinement. 

To summarise: We presented an analysis of quenched QCD in the co- 
variant Landau gauge. Our results show that infrared slavery is at work, 
though in a different fashion than has been assumed in many previous inves- 
tigations. It is not a gluon propagator with a l/p 4 -behaviour in the infrared 
that confines quarks. Instead the gluon propagator is even vanishing at zero 
momentum. However, there is enough infrared strength in the quark-gluon 
vertex to compensate for this: We found a linear rising potential from the 
four-quark Green's function. This potential is triggered by scalar Dirac am- 
plitudes in the quark propagator and quark-gluon vertex. However, since 
A = @2i the potential also contains vector contributions. An answer to the 
old question of scalar vs. vector confinement is therefore a nontrivial dynam- 
ical issue. Corresponding numerical results will be published elsewhere [24]. 
Finally, we wish to emphasise again that if chiral symmetry is restored the 
confining solution disappears. As a result we have uncovered a novel link 
between dynamical chiral symmetry breaking and confinement. 
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